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ABSTRACT 

Recent observations of Type la SNe at redshifts < z < 1 reveal a progressive dimming which has 
been interpreted as evidence for a cosmological constant of J7a ~ 0.7. An alternative explanation of 
the SN results is an open universe with = and the presence of <; 0.1 ^tm dust grains with a mass 
density of fidust ~ (few) x 10 -5 in the intergalactic (IG) medium. The same dust that dims the SNe 
absorbs the cosmic UV/optical background radiation around ~ 1 /im, and re-emits it at far infrared 
(FIR) wavelengths. Here we compare the FIR emission from IG dust with observations of the cosmic 
microwave (CMB) and cosmic far infrared backgrounds (FIRB) by the DIRBE/FIRAS instruments. 
We find that the emission would not lead to measurable distortion to the CMB, but would represent 
a substantial fraction ( > 75%) of the measured value of the FIRB in the 300 — 1000 (im range. This 
contribution would be marginally consistent with the present unresolved fraction of the observed FIRB 
in an open universe. However, we find that IG dust probably could not reconcile the standard fl = 1 
CDM model with the SN observations, even if the necessary quantity of dust existed. Future observations 
able to reliably resolve the FIRB to a flux limit of ~ 0.5 mJy, along with a more precise measure of the 
coarse-grained FIRB, will provide a definitive test of the IG dust hypothesis in all cosmologies. 
Subject headings: cosmology: theory - cosmology: observation - cosmology: Far Infrared Background - 
cosmic microwave background - galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

One of the most remarkable results in recent observa- 
tional cosmology is the detection of type la SNe at cos- 
mological distances. Comparing the Hubble diagram of 
these ~ 50 SNe with the predictions of classical cosmolog- 
ical models strongly favors models with acceleration in the 
cosmic expansion, presumably due to a significant cosmo- 
logical constant (Perlmutter et al. 1999; Riess et al. 1998). 
Given the far-reaching implications of this result, it is im- 
portant to assess the possible systematic effects that could 
mimic the behavior of cosmic acceleration, and thus allow 
cosmologies without a cosmological constant. 

Intergalactic (IG) dust could provide such a systematic 
effect. The amount of dust required to fully account for 
the systematic dimming of SNe at high rcdshift in an open 
universe is fWst ~ (few) x 10~ 5 (Aguirre 1999a; Aguirre 
1999b [A99]). If this quantity of dust had the same com- 
position as interstellar dust in the Milky Way, it would 
cause a significant reddening in the SN spectra, which is 
not observed. However, A99 showed that the removal of 
small dust grains during their ejection from galaxies into 
the IG medium could bias the composite dust cross section 
to cause little reddening for a given amount of extinction. 
A99 presented specific models for the evolution of the IG 
dust that are consistent with the observed lack of redden- 
ing, but provide the amount of extinction to fully account 
for the SNe results in an open universe. 



In the present paper, we study the observational con- 
sequences of the IG dust hypothesis in far infrared wave- 
length bands. As is well known from previous studies, 
IG dust would efficiently absorb the UV/optical flux of 
galaxies and quasars, and re-radiate this energy in the far 
infrared (e.g. Wright 1981; Bond, Carr & Hogan 1991; 
Loeb & Haiman 1997). Although studied in great detail 
in these works, this effect is at present of renewed inter- 
est, owing to recent direct estimates of the global average 
comoving star formation rate (SFR) in the universe (e.g. 
Madau 1999 [M99]), and determinations of the extragalac- 
tic UV/optical background (UVB, Bernstein 1997), as well 
as of the cosmic far infrared background (FIRB; see Fixsen 
et al. 1998; Puget et al 1996). We thus have, for the first 
time, a measurement of both the UV reservoir available 
for dust absorption, as well as the background flux at the 
wavelengths where this energy would be emitted by the 
dust. 

In this paper, we compute contributions to the FIRB 
and CMB arising from IG dust, and compare it to the 
measured values. Our aim is to assess whether the amount 
and type of dust that would explain the SNe results in an 
Q\ = universe is compatible with these new observa- 
tions. 

At the present time, the source of the observed FIRB 
remains unclear, although rapid progress is being made us- 
ing deep observations with SCUBA. Barger et. al (1999a) 
and Hughes et al. (1998) have detected a population of 
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discrete IR sources whose cumulative flux down in flux 
to ~ 2 mJy accounts for at least 20-30% of the FIRB at 
850 /um. The gravitationally lensed sample of Blain et al. 
(1999a, b) goes further down to ~ 1 mJy, but is more un- 
certain because of the lensing model and small number of 
detections. It it therefore currently uncertain whether the 
counts continue to sufficiently faint flux to account for all 
of the FIRB. 

Furthermore, we lack information on the spectra of these 
sources at wavelengths other than 850 /im. Although the 
full FIRB can be explained theoretically in semi-analytic 
galaxy formation models (Guiderdoni et al. 1998), this 
requires the somewhat ad-hoc postulate that the ratio of 
the number of ULIGs (a population of "ultraluminous in- 
frared galaxies") to that of normal optical galaxies in- 
creases rapidly with rcdshift. In summary, according to 
our present knowledge, a fraction 70-80% of the FIRB 
could still be contributed by diffuse emission from IG dust 
- but this may change very quickly. 

As in the case of the FIRB, the nature of the observed 
UVB is unclear. The directly observed values at 0.55 
and 0.8 /zm are 20±10 nWrn^'sr" 1 (Bernstein 1997). By 
comparison, a recent compilation of ground based galaxy 
counts, and a survey of galaxies in the Hubble Deep 
Field (HDF) has yielded an estimate for the evolution 
of the global star formation rate in the universe between 
< z ^ 5. The UVB that follows from a census of these 
galaxies is around ~ 12nWm~ 2 sr _1 (Madau 1999 [M99]; 
Pozzetti et al. 1998), a factor of ~ 2 smaller than the 
observed value. A natural explanation for the difference 
would be if the remaining ~ 50% of the UVB were con- 
tributed by faint, undetected galaxies; however this ap- 
pears to violate limits from the fluctuations in the UVB 
measured in the HDF (Vogeley 1997). Despite these un- 
certainties, the approximate amplitude of the UVB allows 
us to compute useful estimates of the contribution of IG 
dust to the FIRB. 

The rest of this paper is organized as follows. In § 2 we 
summarize our model for the amount and type of IG dust; 
in § 3 we describe our assumptions concerning the redshift 
evolution of both the UV emission and dust production, 
and outline our calculation methods; and in § 4 we discuss 
limits from existing FIRAS observations of the CMB. Sec- 
tion § 5 contains our main results on the contribution to 
the FIRB from IG dust, as well as hopes for future tests 
of the dust hypothesis for the SN results; and, in § 6 we 
summarize our conclusions. Unless otherwise stated, in 
this paper we adopt an open cosmology with total matter 
density ft = 0.2 and Hubble constant h 50 = 1. 

2. IG DUST MODEL 

The model adopted here for IG dust is discussed at 
length in A99; here we provide only a brief summary. The 
model is based on the following method of estimating the 
total IG dust density f2dust(z): we first estimate the to- 
tal metal density Qz(z), then multiply this figure by the 
fraction Fi of these metals which lie outside galaxies, and 
then by the fraction d m of IG metals is contained in IG 
dust. 

As shown in A99, both a direct integration of the SFR 
(with an assumption of the metals produced for each star 

lr The effects are, of course, somewhat dependent upon the grain i 



formed) and fossil evidence from clusters indicate that 
fl z (z £ 0.5) « (2.5 - 5) x 10~ 4 (adjusted from A99 for 
h 50 = 1). These estimates follow from conservative as- 
sumptions: that stars in cluster galaxies have the same 
IMF as in field galaxies, that there is not a dominant pop- 
ulation of unobserved galaxies, etc.; see A99. 

Measurements of the metallicity of intracluster gas indi- 
cate that a fraction Fj w 75% of a typical cluster's metals 
lie in the intracluster medium (e.g. Renzini 1997), pre- 
sumably removed from the galaxies by some combination 
of winds, dust expulsion, ram-pressure stripping, and tidal 
disruption/merging of galaxies. While ram-pressure strip- 
ping and mergers may be more effective at removing met- 
als from cluster galaxies, it is nevertheless likely that the 
figure Fj « 75% applies to the field galaxies, since a much 
smaller value, together with the estimate of Viz <; 10~ 4 , 
would imply that field galaxies have mean metallicity sev- 
eral times solar - contrary to observation. 

To estimate the fraction d m of IG metals in dust, we 
assume a value of ~ 0.5 for metals leaving the galaxy, as 
applies to the interstellar medium of both typical galax- 
ies and perhaps even to damped Ly-a systems (Pei, Fall 
& Hauser 1998). Some fraction (1 — /) by mass of this 
dust must be destroyed either during ejection or in the IG 
medium, giving d m « 0.5/. Combining these figures, we 
estimate an IG dust density at z <; 0.5 of 

fidust* (9.4- 18.8)/ xl0~ 5 . (1) 

To fix both / and the dust properties, we shall adopt the 
two-component Draine & Lee (1984) [DL] dust model 1 
of silicate and graphite spheres with size distribution 
N(a)da oc a" 3 5 , 0.005 /im < a < 0.25 ^m (Mathis, Rumpl 
& Nordsieck 1977 [MRN]) as representative of galactic 
dust. We further assume that IG dust differs from MRN 
dust due to the removal - by the selective ejection of large 
grains and/or the selective destruction of small grains by 
sputtering - of the small-size end of the grain-size distribu- 
tion. In this paper we use absorption and scattering curves 
as calculated in Laor & Draine (1993), and the truncated 
MRN distribution with a m ; n = 0.1 /im; The minimal grain 
size (which corresponds to / w 0.4) is chosen to give dust 
which is grey enough not to over-redden the supernovae 
(A99). Note that our results should not be sensitive to our 
adoption of the Milky Way dust cross section - e.g. LMC 
dust differs from Milky Way dust mainly at short wave- 
lengths, where the opacity has anyway been modified by 
removing the small grains. Silicate and graphite grains are 
assumed to have equal mass densities, and the total dust 
density is chosen, according to the cosmology, as that suffi- 
cient to reconcile the chosen cosmology with the supernova 
results. For example, = 0.2 requires ~ 0.15 — 0.2 mag of 
extinction at z <~ 0.5 and hence f^dust ~ (5.5 — 7) x 10~ 5 
to be consistent with the observations. 

3. METHOD OF CALCULATION 

In estimating the contribution from IG dust to the FIRB 
and the CMB, we follow the methods of Wright (1981) 
(see also Loeb & Haiman 1997). We define the comov- 
ing number density of photons at redshift z and comoving 

; see A99 for discussion. 
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frequency v as 

N v {z) 



-in 



hc(l + zf 



(2) 



where 3 v (z) is the usual specific intensity of the back- 
ground radiation field in physical (non-comoving) units 



of erg cm 



2 s- 1 



Hz sr . The evolution of N v {z) with 



redshift is then described by the cosmological radiative 
transfer equation 



dN v (z) cdt 



dz 



dz 



\jtot,v(z) - a v (z)N v (z)\ 



(3) 



where jtot,v — 3±,v + 3&ust,v is the total comoving emission 
coefficient, including both the direct light of UV sources 
(see equation H below), as well as thermal dust emis- 
sion jdust,^ from dust. As discussed in section § 2 above, 
we assume that dust is composed of graphite and sili- 
cate grains, so that the total dust absorption coefficient, 
a v = asi,v + ctGr,i/ is given in units of cm -1 by 

a v (z) = PSi(2)«Si,i/(l+z) + PGr(z)K Gr ^ {1+z) (4) 

for a total mass density of psi+Pa- = (^Si+^Gr)Pcrit,o(l + 
z) 3 in dust, where the re's are the dust opacities in units of 
cm 2 g _1 . Note that by assumption, figi = ^Gr- The total 
comoving dust emission coefficient is then 



jdust.K = jsi,v + 3Gx,v 
8tw 3 , 



asi,u(z) 



c 3 L exp[/u/(l + z)/k B T Si ] -1 

OLGt,v{z) 



exp[^(l + z)/k B T Gl ] - 1 



}■ 



(5) 



(6) 



The silicate and graphite grains are assumed to be sep- 
arately in thermal equilibrium with the total background 
radiation (CMB + UV sources), so that the grain tem- 
peratures Ti are obtained by equating the total amount of 
absorption and emission: 



(7) 



dv a- hV N v (z) = / dv j-^ 



where i=Si or Gr. 

The infrared background spectrum arises from the re- 
emission of UV light absorbed by dust. The two main 
factors determining the F1RB are therefore the evolution 
of the UV emissivity j+ .„ from stellar sources (determin- 
ing the dust temperature), as well as the evolution of the 
dust density, ^Wst = ^Si + ^Gr - determining the overall 
amplitude of the emissivity. For the UV production rate, 
here we adopt the global average star-formation rate as 
determined recently by Madau et al. (1998). We ignore 
any additional UV flux from quasars: at the redshifts of 
interest here (z <> 5), the contribution to the UV back- 
ground from the known population of quasars is less than 
20% (M99). Note, however, that a population of faint, un- 
detected quasars could contribute substantially to the UV 
background (Haiman & Loeb 1998). 

For simplicity, we assume further that the average spec- 
trum from the stellar UV sources is described as that of 
black-body radiation at a temperature of T = 9000K. This 



is a reasonably good approximation to the continuum in 
the composite spectra found in the population synthesis 
models of Bruzual & Chariot (1996), although the pure 
black-body spectrum is somewhat narrower. The black- 
body also falls below the UV background flux limit at the 
shortest wavelengths (A ~ 0.1 /im); see also § 3. Our re- 
sults below depend primarily on the total amount of stel- 
lar light absorbed by the dust, and are insensitive to the 
precise shape of the UV spectrum (see detailed discussion 
below) . With the above assumptions, the stellar emissivity 
is given by 

h, u (z) = A uv p+(z)B H1+z) (T = 9000K), (8) 

where B v is the Planck function, and p*{z) is the 
star formation rate in Madau et al. (1998), converted 
from their standard fl = 1 cosmology to the open Q = 
0.2 model assumed here, using the multiplicative factor 
(dVo/df )/ (dv\/df A ), where d v and di are the volume el- 
ements and luminosity distances in the two cosmologies 
(Eisenstein 1997). 

The overall normalization constant A uv is somewhat un- 
certain. M99 argues that the total UV produced by stars 
is a factor ~ 2 smaller than the recent measurement of the 
extragalactic UV background, vJ v = 20 nW vaT 2 sr -1 at 
0.55 /im by Bernstein (1997). It is not clear what the rea- 
son for this discrepancy is; the shortage of UV could be 
provided by additional faint galaxies or quasars (however, 
see Vogeley 1997). Here we take the simplest approach, 
and adjust the normalization constant A uv so that the 
measured value of 20 nW raT 2 sr -1 for the UV background 
at 0.55 /xm is exactly reproduced (taking into account dust 
absorption). Given the shortage of UV provided by the 
known stellar population, we will also examine below a 
model normalized to 10 nW m -2 sr -1 , the lower end of 
the la range quoted by Bernstein (1997). This also lies 
above the lower limits on the UV background derived from 
galaxy counts in the Hubble Deep Field (Pozzetti et al. 
1998; see figure ||). Note that we have assumed a specific 
spectral shape and redshift-evolution of the UV emissiv- 
ity. We again emphasize that our results are insensitive 
to the detailed shape of the UV flux around 1 /xm, and 
depend only on the total amount of UV light absorbed; 
this is demonstrated by a model we examine below with 
a more realistic spectrum obtained from population syn- 
thesis models (Bruzual & Chariot 1996) and models with 
SF histories different from M99. As mentioned above, the 
difference between the background inferred from the ob- 
served star-formation rate (i.e. by direct integration of /)+ 
and an empirical UV flux per stellar mass ratio), and the 
measured value of 20 nW m~ 2 sr -1 could be provided by 
undetected quasars, or other yet undiscovered sources. As 
long as the redshift evolution of these sources does not de- 
part significantly from the range of our assumed p^(z)'s, 
our results below would remain valid. 

Next, we assume that the same stellar population that 
produces the UV background is the source of IG dust. 
Under this simple assumption, the mass density of dust 
increases in proportion to the mass density in stars: 



OdustO*) = Adust / dtp^z), 

Jt(z = W) 



(9) 



and we assume further that f^si = ^Gr = ^dust/2 (as dis- 
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cussed in § 5 the exact value of this ratio is not crucial). 
The normalization constant ^dust is chosen by requiring 
the total optical depth due to dust absorption + scatter- 
ing between z — and z — 0.5 to be ~ 0.15. As we 
argued in section § 1 above, this is the amount needed 
to account for the SN results (A99). This normalization 
corresponds to Odustt 2 = 0) ~ 5 x 10~ 5 . Adding to the ar- 
guments given in § 1, we note that this dust density fits in 
comfortably with the amount of dust expected to be pro- 
duced in normal stars. For instance, assuming that each 
type II SN (M* > 8M ) yields 0.3M Q of dust, adopting a 
Scalo IMF, and normalizing f2*(z = 0) = 0.004; we obtain 
f^dustO = 0) = 5 x 10~ 5 (cf. Loeb & Haiman 1997). This 
estimate does not account for dust trapped in galaxies or 
destroyed, but is also a lower limit since it only includes 
the dust formed in supernovae, and not that formed in (for 
example) dense clouds or carbon stars. 

We note that although stellar dust production cannot 
occur instantaneously with star-formation, the relevant 
time-scales are likely to be much shorter than the Hub- 
ble time at redshifts of interest. Type II SNe produce 
dust in a few million years; while the lifetime of carbon 
stars (believed to be dominating the dust budget in the 
Milky Way, see Gehrz 1989) is still only -1 Gyr. The 
ejection timescale for dust leaving spiral galaxies is proba- 
bly between 100 and 1000 Myr (A99). We have found that 
incorporating a delay of ~1 Gyr between star formation 
and dust production into equation ^| does not change our 
results below by more than a few percent. 

Equations together with the initial condition that 
N v at some high redshift is equal to the pure CMB with 
T cm h = 2.728(1 + z)K, determine the redshift evolution of 
the full background spectrum. 

4. EXISTING FIRAS LIMITS 

Any existing IG dust component will process some 
UV/optical energy into the FIR/microwave. As a result, 
some (perhaps vanishingly small) fraction of the observed 
CMB is due to emission from IG dust. This contribu- 
tion has been studied a number of times previously. Some 
authors have investigated the possibility that the CMB 
is entirely dust-processed emission from astrophysical ob- 
jects; see, e.g., Layzer & Hively (1973), Wright (1982), 
Wollman (1992) and Aguirre (1999c) on 'cold big-bang' 
models, or Wickramasinghe et al. (1975), and Hoyle & 
Wickramasinghe (1988) on steady-state models. In cold 
big-bang models, thermalization takes place at very high 
z after energy emission by Population III objects. Both 
cosmologies require grains of a type that can maintain a 
temperature very close to that of the CMB. 

If the grains have a temperature only slightly different 
from that of the CMB, they will be too cool to contribute 
to the FIRB; however, spectral distortions of an initially 
blackbody CMB will occur (e.g. Rowan-Robinson, Negro- 
ponte & Silk 1979; Hawkins & Wright 1988; Bond et al. 
1991). More recently, this effect has been used by Loeb & 
Haiman (1997) and by Ferrara et al. (1999) to limit the 
density of IG dust. We shall perform a similar analysis 
using the method and assumptions outlined in § ||, and 

2 In light of the uncertainty in the FIRB, it seems that the quoted 
must be used with caution. Rather large spectral distortions can be 
spectral shape. Assuming conservatively that the FIRB as detected by 
to a ^-parameter in the full 60-600 GHz interval of 2 X 10 -4 , over an 



with the dust model described in § 1. 

Models which distort the CMB at long (3> 1mm) wave- 
lengths can be easily constrained since there is a clear 
model for the CMB (perfect blackbody). But in the ab- 
sence of an equally unique and well-defined model for the 
FIRB, it is much less clear how to limit FIR emission that 
does not actually exceed the observed FIR background. 
The procedure adopted here is as follows: following the 
methodology of Fixsen et al. (1996), we fit the results 
from 500 — 5000 fxm by a blackbody -B„(T cm b) plus a 'uni- 
form dust component' parameterized by an emissivity and 
temperature in a v 2 B^{T) spectrum. These fitted compo- 
nents are shown in dotted lines in Figure ^. The spectrum 
with the CMB only subtracted is shown as a solid line; 
this makes clear the contribution of IG dust emission to 
the FIRB. The residuals with both components subtracted 
are shown in the second panel, along with the FIRAS fre- 
quency range (vertical dotted lines) and the limit on r.m.s. 
deviations from blackbody (horizontal dotted line). 

An interesting note about this general procedure is the 
following. Our calculations show that the fitted CMB tem- 
perature may be close to, but not exactly the same as the 
actual initial CMB, and this small discrepancy affects the 
measured slope of the FIR emission at A ;> 850 /xm, since 
the FIRB is < 1/1000 of the CMB flux there. For ex- 
ample, fitting the CMB plus v 2 B v gives a fitted CMB 
temperature of 2.72800K, exactly as it should be. But 
if we fit a CMB plus i/°- 64 B„ (the shape of the FIRAS 
fit), we find T cmb = 2.72795K. The CMB is then slightly 
under-subtracted, and the residual FIRB closely follows a 
/„ oc j/ 3 - 64 shape at long-A. This suggests that accurately 
determining the slope of the FIRB long-wavelength tail 
requires knowledge of T cm b to much higher accuracy than 
available from COBE data. For this reason, we find it is 
not useful to compare the FIRB slope predicted by models 
to the FIRAS fit, for A £ 850 fim. 

Table gives the r.m.s. deviation from blackbody for 
the computed spectra, with the mean taken over the 
500 /im < A <J 5000 fim wavelength range analyzed by 
Fixsen et al. (1996) and divided by the peak of B v {T cm h). 
The limit given by Fixsen et al. 1996 on this number is 
5 x 10 -5 . The result of this analysis is that the dust emis- 
sion would not cause detectable distortions to the CMB for 
the dust model described, despite the fact that the amount 
of emission appears to violate the quoted limits on the y 
parameter. 2 The residuals are, of course, even smaller if 
the 'uniform dust component' is fit using a free index in 
the power law. 

On the other hand, modification of the long-wavelength 
opacity of the dust can be important, since dust types with 
higher FIR opacity have lower equilibrium temperatures 
which may lead to excessive emission into the FIRB/CMB. 
Since the dust emission tends to be well described by the 
fitted v 2 B v isotropic dust component, this emission does 
not tend to leave excessive distortion in the CMB residuals. 
Rather, it leads to a fit of the isotropic dust component 
incompatible with that found by FIRAS. Accordingly, we 
discuss these models in the next section, which treats the 



limits on r.m.s. distortions and the {/-parameter from the COBE group 
'hidden' in the microwave tail of the FIRB which has a poorly defined 
FIRAS exists in the 150-600 GHz range (and zero outside), it corresponds 
order of magnitude above the quoted upper limit. 
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dust contribution to the FIRB and the resulting limits on 
the dust models. 

5. IG DUST CONTRIBUTION TO THE FIRB 

In this section, we analyze the contribution by IG d ust 
emission to the FIRB in a variety of models. Section § [Tl] 
presents and discusses our fiducial model. Section § 3 lists 
possible variations in the parameters, and discusses the 
effects of these variations, and section § 4 discusses how 
resolution of the FIRB into discrete sources by future ex- 
periments would improve the constraints. 

5.1. Fiducial Model 

In our fiducial model, the cosmology is open, with 
Q = 0.2 and = 1. The star formation rate is taken 
from M99 and shown in A99, Figure 1. The dust has 
equal parts by mass in silicate and graphite grains with 
Ornin = 0.1 /im, and density proportional to the integrated 
SFR, normalized to f2dust(z = 0) = 5.4 x 10 -5 , which gives 
0.15 mag of extinction to z = 0.5 at (observed) 0.66 /mi. 
The galaxy spectrum is a 9000 K blackbody, normalized 
to give 20nWm~ 2 sr~ 1 after dust processing at 0.55 /mi. 

The predictions of the fiducial model appear in Figure |^ 
and Table Q. The fi gure shows the complete final spectrum 
with the CMB (fit using the method of § 2) subtracted. For 
comparison we also include the UVB lower limits (Pozzetti 
et al. 1998) and detections (Bernstein 1997), the DIRBE 
FIRB upper limits and detections (Hauser et al. 1998), 
and the FIRAS detections (Fixsen et al. 1998) with ±la 
uncertainties. It is apparent that the dust emission con- 
tributes significantly to the long-wavelength flux - in fact 
it can account for all of the 850 /im flux given by the FIRAS 
measurements - but fails to reproduce the shape or ampli- 
tude of the entire FIRB. This is specified quantitatively in 
Table |l|, where the fifth and sixth columns list the numer- 
ical fraction of the FIRAS FIRB measurements that the 
dust emission would contribute at 200 and 850 /.tm. These 
numbers should be multiplied by 0.76 to yield the fraction 
of the la upper limit on the FIRB, or 1.25 to account for 
the minimum 20% resolved fraction. Column seven lists 
the y-parameter represented by the FIR emission, defined 
by y = 0.25(wtot/u cm b - 1), where u to t and M cm b are the 
total energy density and the energy density of the CMB 
alone, both evaluated in the 60-600 GHz frequency range. 
This gives a measure of the total dust energy output. The 
last column gives the r.m.s. residual after both CMB and 
the v 2 B l/ fit to the FIR emission are subtracted, as dis- 
cussed in § 2. 

The next section analyzes the possible variations upon 
these predictions resulting from changes in the model. 
Some of the variations can be ruled out because they pre- 
dict a FIRB contribution above the FIRAS limit; these 
models could be saved only if the FIRAS team has over- 
subtracted the foreground contribution. Stronger limits 
can be obtained on the various models by considering their 
contribution to the unresolved fraction of the FIRB, af- 
ter the contributions by known discrete sources are sub- 
tracted. These constraints will be discussed in section 4 
below. 



5.2. Variations on the model 

Several variations were performed to test the robust- 
ness our model assumptions. First, raising Hq from 50 to 
80kms _1 mpc _1 (model H80) made very little difference. 
Second, an SFR which is flat for z > 3, but unchanged 
for z < 3 was employed (model FSFR). This choice is 
motivated by the "effective" global average SFR (M99) 
required to reionize the universe by z — 5. This increases 
the dust contribution to the final spectrum by about 35%, 
and the 850 /im contribution by 45%. The change occurs 
because relatively more of the UV energy is emitted at 
high z, where the dust is more dense, and emits into longer 
(rest-frame) wavelengths. 

Next, we tested the adequacy of our assumption of the 
single-temperature blackbody galaxy spectrum, by using 
a more realistic spectral template. We utilize the Bruzual- 
Charlot (1996) starburst population synthesis model with 
a Scalo IMF to compute the shape of spectrum as a func- 
tion of the age of the population. This time-dependent 
template is convolved with the SFR from M99 to find the 
rcdshift-cvolution of the frequency dependent UV emissiv- 
ity. The full spectrum under these assumptions is shown 
by the curve labeled "BC" in Figure [|. The peak in the 
UVB from the direct galactic emission is broader than in 
our fiducial model, and extends to higher energies. As 
mentioned above, this procedure results in a UVB at z = 
that falls short of the measured value by a factor of ~ 2. 
Overall, the contribution to the FIRB is not effected sig- 
nificantly: it is reduced only by ~ 10% relative to our 
fiducial model at 850 /im (cf. Table |). 

It must be noted that the above procedure is not strictly 
correct, as we have not corrected the galactic templates 
for reddening by dust internal to the galaxy. This red- 
dening could be important in the present context, since 
it makes the UV peak steeper in the 0.1 — 1 /im range, 
and allow a broader UVB in this wavelength range that 
reaches ~ 20nWm~ 2 sr~ 1 at 0.55 /mi, while extending to 
energies higher than assumed in our fiducial model. Here 
we do not attempt to derive an accurate reddening correc- 
tion (see instead Madau et al. 1998). Rather, our aim is 
to assess the maximum UV flux that could, in principle, 
be present in the 0.1 — 1 /im range, that is consistent with 
both the Voyager upper limit (Murthy, Henry & Holberg 
1998) and the 0.55 /im detection. 

To compute this maximal UV flux, we reddened the 
galactic template spectrum, assuming the cross section of 
Milky Way dust (DL), with a total (unrealistically high) 
maximum optical depth of ~ 4. We then re-scaled the 
UV efficiency by a factor of ~10 to retain an intensity of 
~ 20nWm~ 2 sr~ 1 at 0.55 /mi. The resulting spectrum is 
shown by the curve labeled "BC2" in Figure g. This spec- 
trum fits in with both the Voyager and Bernstein (1997) 
data points, and reasonably represents the maximally al- 
lowed UV emission. It is interesting to note that a simul- 
taneous fit to these two data-points requires a very steep 
spectrum (20nWm~ 2 sr _1 at 0.55 /mi to 0.6nWm~ 2 sr 
at 0.1 /xm, implying a slope of I ~ A 2 ) , which can only 
be derived from the Bruzual-Charlot template by postu- 
lating an exceedingly high reddening. One resolution of 
this paradox might be that the true value of the UVB is 
closer to 10nWm~ 2 sr _1 (as argued by Vogeley 1997 in a 
different context). We find that with this maximal emis- 
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Table 1 

intergalactic dust emission models 



Model 


V ClllClLlWll 


Trv 
± Or 


J Ol 


f 

J200 


f 

J85Q 


7/T7TTD T3 / 1 

</r 1KB / ±KJ 


rms/10~ 6 


Fid. 


- 


9.9 


7.3 


0.14 


1.02 


2.0 


3.7 


H80 


Hq = 80kms _1 mpc -1 


9.9 


7.3 


0.15 


1.10 


2.1 


4.0 


FSFR 


SFR flat for z > 3 


9.9 


7.3 


0.16 


1.45 


2.7 


7.4 


BC 


Galaxy spectra from Bruzual-Charlot 


9.3 


7.8 


0.08 


0.89 


1.7 


2.9 


BC2 


Bruzual-Charlot + extreme reddening 


10.5 


7.9 


0.24 


1.26 


2.5 


4.6 


MG 


a = 0.1 /jm grains, Qd = 4.5 x 1CP 5 


10.2 


7.1 


0.17 


0.96 


1.9 


3.4 


LG 


a = 0.25 fim grains, fid — 7.6 x 10~ 5 


9.2 


7.4 


0.08 


0.85 


1.6 


3.1 


LUV 


UVB of lOnWm^sr- 1 at 0.55 /un 


8.8 


6.5 


0.05 


0.72 


1.3 


2.5 


SCDM 


n = i, n d (z = o) = 1.25 x io~ 4 


9.9 


7.3 


0.30 


1.94 


3.8 


6.7 


SCDMb 


SCDM, but UVB of lOnWm^sr 1 


8.8 


6.5 


0.10 


1.33 


2.5 


4.4 


G05 


k flat to Ao = 0.5 /im 


17.3 




0.86 


0.45 


1.1 


1.1 


G2 


K flat to Ao = 2 /im 


12.8 




0.88 


1.56 


3.5 


4.6 


G5 


K flat to Ao = 5 fim 


9.7 




0.53 


3.97 


7.8 


12 


G10 


k Hat to Ao = 10 fim 


7.7 




0.24 


7.71 


13.4 


19 


L1.0 


k(A > 100 fim) oc A" 1 


7.9 


5.5 


0.038 


2.59 


4.4 


18 


LI. 5 


k(A > 100 fim) oc A -15 


9.0 


6.4 


0.083 


1.71 


3.1 


9.5 


L2.5 


k(X > 100 fim) oc A -2 ' 5 


10.6 


8.0 


0.19 


0.65 


1.3 


1.7 


GRA 


graphite, <z m i n = 0.05 /im 


9.8 




0.21 


1.29 


2.6 


4.4 


S2 


1:2 graphite: silicate by mass 


9.9 


7.3 


0.12 


1.00 


1.9 


3.7 



sion, the FIRB contribution increases somewhat, but the 
change is once again not significant, only ~ 40% relative 
to the galactic model without reddening, over-predicting 
the 850 fim FIRB by ~ 26% (cf. Table]]]), but still well 
within the +ler uncertainty. 

Third, we have tested our assumption that we may treat 
the grain-size distribution as a single component, by con- 
sidering grains of a single grain size a, for a = 0.1 /im 
(model MG) and a = 0.25 /im (model LG). These ex- 
tremes change the z — temperature by 0.1— 0.7K and the 
850 /im flux by < 4—15%, demonstrating that the approx- 
imation is not very important. The test also shows that 
the results are not sensitive to the grain size distribution, 
as long as all of the grains are relatively large (a ^ 0.1 /mi). 
The exact graphite:silicate ratio is also fairly unimportant; 
models GRA and S2 have graphite:silicate ratios of 1:0 and 
1:2 respectively. The former corresponds to a scenario in 
which silicates are preferentially destroyed (see A99); the 
latter is an extreme of reasonable graphite:silicate ratios 
for two-component model, as shown in the recent compu- 
tations of Weingartner & Draine 1999. Neither is drasti- 
cally different from the fiducial model, and the S2 model 
is almost identical. 

Finally, we check our assumption that the dust temper- 
ature follows from equilibrium with a homogeneous back- 
ground radiation bath, because dust near concentrations 
of radiation could have somewhat higher temperature. To 
determine whether this might be important, we have com- 
puted the critical distance r cr it(z) (from center of a model 
galaxy) at which a dust grain would feel equal contribu- 
tions from the cosmic UVB intensity vj v {z) and from the 



nearby galaxy, i.e. 

in(l + z) 3 vJ„ = -^-, (10) 

where the galaxy luminosity L ga \ ~ 10 46 ergs _1 (corre- 
sponding to a 10 12 M© starburst galaxy, with the spectral 
model of Bruzual & Chariot 1996) is constant and vJ u (z) 
is proportional to the integrated comoving SFR and has 
the value 20nWm~ 2 sr~ 1 at z = 0. The result is that 
r crit <5 60kpc for all z and r cr j t <J 30kpc for z > 1. In 
A99, it was argued that IG dust must be ;> 70kpc from 
its progenitor galaxy to be uniform enough to be in accord 
with limits on the dispersion in supernova brightnesses. 
Thus even near extremely bright galaxies, dust of the A99 
model should have temperature dominated by the isotropic 
background rather than nearby galaxies. 

As discussed in there are still sizeable uncertainties in 
the normalization of the UV spectrum. The la lower limits 
from Bernstein (1997) roughly coincide with the lower lim- 
its given by integrated counts (M99) at ~ 10 nW nr 2 sr _1 
at 0.55 /im. We have employed this normalization (model 
LUV), and find that this lowers the 850 /im flux by about 
30% (the decrease in energy is somewhat offset by the 
lower temperature of the dust). 

Although most models considered predict an 850 /im 
flux comparable to that detected by FIRAS, some models 
exceed the FIRAS upper limit. To reconcile the supernova 
measurements with a closed, matter-dominated cosmology 
- such as the standard cold dark matter model (SCDM) - 
an extinction of Ay (z ~ 0.5) ~ 0.35 mag is required (A99), 
in turn requiring fld{z = 0) ~ 1.25 x 10 -4 for h$a = 1. In 
the present model this predicts an 850 /im flux which ex- 
ceeds the 1 — o bound of the FIRAS detection (see the 
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SCDM model in Table Q.) With the lower UV normaliza- 
tion, however, the emission is marginally compatible with 
the errors, (model SCDMb in Table Q). If the higher UV 
normalization holds, these results make it very unlikely 
that IG dust could reconcile SCDM with the supernova 
observations, even if the large necessary quantity of dust 
existed (cf. A99). 

The present calculations can also constrain the proper- 
ties of intergalactic dust if it is assumed to be responsible 
for the supernova dimming. Very grey dust, for instance, 
absorbs more UV/optical flux for a given V-band extinc- 
tion (and can hence overproduce the FIRB). Both very 
grey dust and 'standard' dust with high IR emissivity can 
also have lower equilibrium temperature, which overpro- 
duces the FIRB at long wavelengths. To limit these sorts 
of dust, we have considered two illustrative types. The 
first is dust that is grey out to some wavelength Ao, then 
falls as A~ 2 : 

1 



K\ OC 



[1 + (A/A ) 2 



(11) 



This is conservative in the sense of ensuring the maxi- 
mum reasonable temperature for a given Ao . Models (G05- 
G10) are of this type, and show that for Ao ;> 3 /im the 
FIRB is overproduced. It is interesting to note, how- 
ever, that Ao ~ 0.5 models actually produce a FIRB fit- 
ting the observations rather well over the entire frequency 
range. Conducting needles have an absorption spectrum 
which can be roughly modeled by eq. [ll| (Wright 1982; 
Wickramasinghe & Wallis 1996; Aguirre 1999a), and using 
Wright's 1982 RC model with a resistivity of ~ 10~ 15 s, we 
can rule out needle models with length-to-diameter ratio 
L/d^6. 3 

The second dust type is 'standard' DL/MRN dust, but 
with long-wavelength (A J> 100 /im) emissivity given by 
k\ cx A _Q with a^2. Fluffy/fractal grains can give a ~ 1 
(e.g. Wright 1987; Stognienko et al. 1995), and a generally 
depends upon the optical properties of the grain material. 
Dust with a < 2 generally has lower temperature, and the 
calculations (see models L1.0, LI. 5 in Table ll) show that 
such dust will overproduce the 850 /jm FIRB (at the la 
level) for a < 1.5. On the other hand, if dust had a > 2, 
model L2.5 shows that the higher resulting temperature 
would lead to less flux at 850 /im. The comoving dust 
temperature Id ust /(1 + z) for the fiducial model and mod- 
els GR05 and L1.0 is shown in Figure |l|. All computed 
models fall within the range spanned by the GR05 and 
L1.0 models. 

In summary, using only the FIRAS limits on the FIRB, 
we can with some confidence rule out SCDM models, 'very 
grey' dust models with Ao > 3 /im, and dust models with 
lower law index < 1.5 in their long- wavelength emissiv- 
ity. The next section discusses stronger constraints which 
derive from the fraction of the measured FIRB which is 
known to be due to discrete sources. 

5.3. Resolving the FIRB: Tests and future observations 

When constraining models using their predicted IG dust 
contribution to the FIRB, we have so far considered the 
full value of the FIRB. However, as mentioned above, 
we already know that the entire FIRB cannot be com- 

3 Needles with high L/d and high conductivity, such as the iron whiskers of Hoyle & Wickramasinghe (1988), can maintain a temperature 
very close to that of the CMB and are not constrained by the present calculations; see also Aguirre (1999c). 



prised of IG dust emission alone. On the one hand, the 
IG dust models generally fail to produce the short wave- 
length fluxes; on the other hand, observations have al- 
ready resolved a non-negligible fraction ( ;> 25%) of the 
FIRB into discrete sources. If 20 - 30% of the FIRB is 
accounted for by resolved sources, the numbers in column 
six of Table |l| should be multiplied by 1.25-1.5, as only 
70 — 80% of the measured FIRB can be IG dust emis- 
sion. To use instead the la upper limit on the detected 
FIRB, but taking into account the resolved fraction, col- 
umn six should be multiplied by 0.95 — 1.14. The cal- 
culations therefore indicate that resolving the FIRB fur- 
ther could provide significantly stronger constraints on the 
models. It must also be noted, however, that there is a 
~ 10 — 20% uncertainty in the SCUBA flux calibration, 
and that the quoted errors on the measured FIRB are re- 
ally rough estimates of systematic uncertainty reflecting 
only differences between the fits resulting from different 
foreground subtraction method. For comparison, Guider- 
doni et al. (1997) give an allowed FIRB intensity range at 

- 850 /im of - 0.2-1.5 xlO^nWm^sr 1 at 1-a. These 
uncertainties correspond to roughly four to five times those 
quoted by Fixsen et al. 

At the present time, a handful (~ 30) of objects with 
fluxes S above S > 2 mJy at 850 /im have been detected by 
SCUBA. These data have been used to derive the empirical 
number counts n(S) of sources, as a function of apparent 
flux. At the limiting flux of S = 2 mJy, the total contribu- 
tion from these objects is a fraction 20 — 30% of the FIRB 
value (Barger et al. 1999a). As emphasized by these au- 
thors, the faint-end counts are very steep {n(S) oc S 1-3 2 ). 
Extrapolating at the faint end of the number counts to 

— 0.5 mJy would then account for the entire measured 
background. These conclusions are further supported by 
the somewhat less reliable detections of faint sources down 
to « 0.75 mJy: the gravitationally lensed sample of Blain 
et al. (1999a, b) extends to this limit, albeit with only a 
four sources fainter than ~ 2 mJy. 

In light of these existing measurements, and especially 
because of the steepness of the faint-end number counts, 
tightening the constraints obtained in this paper requires 
relatively small further improvements in the observations. 
In particular, with a faint-end slope of -3.2, the flux limit 
needs to be lowered only by a factor of four, to ~ 0.5 mJy, 
to detect the sources that account for the whole of the 
FIRB, provided these sources exist. If the faint sources do 
not exist, then these measurements will find the turn-over 
and flattening in n(S), and yield the actual fraction of 
the FIRB not accounted for by these sources. If the faint 
sources are detected, the uncertainty in their total cumu- 
lative flux would be relatively small, since these sources 
are expected to be numerous. 

A complication we have ignored in this analysis is the 
following: most IG dust emission occurs at 1 <; z <; 3, 
where a tentative redshift analysis (Barger et al. 1999b 
and references therein) indicates that most of the SCUBA 
850 fim sources lie. The source fluxes are derived using 
a 30 arcsec aperture, which at z > 1 corresponds to a 
physical size of ;> 250/i 65 1 kpc. We have assumed for sim- 
plicity that the iC dust is effectively uniform, but at z > 1 
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significant quantities could be within 125 kpc of galaxies, 
so there is a danger that the SCUBA derived fluxes in- 
clude some emission from what we have termed IG dust. 
Upcoming sub-mm surveys with better angular resolution 
should be able to clarify this ambiguity. 

In summary, with a factor of four improvement over the 
current SCUBA detection threshold, we would either see 
explicitly the need for something else to explain the FIRB, 
or else a fairly tight constraint on the IG dust will be pos- 
sible. Forthcoming instruments such Bolocam and SIRTF 
will reach the required sensitivities at 160 fim and 1100 /xm 
(see Table 1 in Blain 1999). The constraints that can be 
obtained then might then be limited by the accuracy of the 
absolute measurement of the coarse-grained FIRB. The 
High Frequency Instrument (HFI) on the future Planck 
satellite 4 , covering the 350 — 3000 /im range at several in- 
termediate frequencies, will be able to greatly improve our 
existing knowledge of the unresolved component. Finally, 
as shown by Haiman & Knox (1999), angular correlations 
in the FIRB are expected to be at the few percent level, 
and depend strongly on the nature of its sources. Future 
instruments such as BLAZE and FIRB A Twill be able to 
measure these correlations, providing another diagnostic 
that distinguishes between IG dust and discrete sources. 

6. CONCLUSIONS 

In this paper, we studied the far infrared emission from 
the type and amount of IG dust necessary to explain 
the recent Hubble diagrams, derived from observations of 
Type la SNe at redshifts < z < 1, in cosmological models 
without a cosmological constant. In particular, we com- 
puted the contribution from the IG dust emission to the 



value of the FIRB recently measured by the COBE satel- 
lite. We investigated a broad range of models, and focused 
on the wavelength of 850 /im, where the largest fraction 
of the FIRB is presently resolved into discrete sources, 
thus yielding the strongest constraints. The current con- 
straints are limited by ~ 10 — 30% errors in each of the 
following: the dust model, the measured FIRB flux, the 
SCUBA counts and the SCUBA calibration. 

Our results show that the IG dust emission is consistent 
with the spectral distortion of the CMB allowed by present 
COBE data, but may contribute nearly all of the unre- 
solved fraction of the FIRB in the 300 — 1000 fim range. 
In a few specific models, this contribution is sufficiently 
large to render those models implausible, including Stan- 
dard CDM and models with very grey dust, or dust with 
unusually high IR emissivity. In the (perhaps most inter- 
esting) case of an open universe with a low matter density 
(fi = 0.2), we find that the contribution is still within 
the experimental uncertainty. Assuming that 20% of the 
FIRB is accounted for by the discrete SCUBA sources, 
the IG dust contributes up to ~ 90% of the +la limit on 
the unresolved fraction. Future observations of the far in- 
frared background by Planck, and its discrete constituents 
by SIRTF and Bolocam will provide a definitive test of the 
IG dust hypothesis. 
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was supported in part by the National Science Foundation 
grant no. PHY-9507695. 
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Fig. 1. — Evolution of dust temperature for silicate (dotted) and graphite (solid) grains in our fiducial model (F), a dust model with 
long-wavelength opacity oc A -1 (L1.0) and a model with opacity that is constant for A < 0.5 (im and falls as A~ 2 for longer A (GR05). The 
dust temperatures in all other models fall in-between the two extreme models shown here. 
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100 1000 10000 

A [fJ.m] 

Fig. 2. — Top: Full z = output spectrum (solid line) for our fiducial model, with the best-fit CMB (higher dotted line) subtracted. 
Also shown are the fitted v 2 B u component (lower dotted line), the FIRAS detected FIRB ±1ct bounds (dashed), DIRBE FIRB detections 
(crosses) and upper limits (down arrows), the UVB detections (triangles) and HDF-derived lower limits (up arrows), and the Voyager FUV 
upper limit (diamond). Bottom: Spectrum with both CMB and v 2 B v component subtracted (solid and dashed lines). Vertical dotted lines 
represent the wavelength range over which the FIRAS distortion limit analysis was performed, the horizontal dotted line indicated the r.m.s. 
of the residuals in the FIRAS data. The features in the 500-5000 fim range indicate the level of numerical noise in our calculations. 
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A [/im] 

Fig. 3. — Output spectra for various models: A model using half the UVB normalization (LUV; dotted), an Q = 1 model (SCDM; dashed), 
and models with the Bruzual-Charlot spectral templates (BC; triple-dot-dashed), and reddened B-C templates (BC2; dot-dashed). The 
fiducial. model is included (solid), as are UVB and FIRB limits as in Figure 0, with FIRAS FIRB detection limits in dark, solid. See text and 
Table M for details on models. 
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X [yiim] 

Fig. 4. — Output spectra for different dust models: A model with long-wavelength opacity oc X~ a for a = 1 (L1.0; dashed) and a = 2.5 
(L2.5; long dashed), and models with dust which is grey out to Ao then falls as A~ 2 , for Ao = 0.5 /an (GR05; triplGidot-dashed) and for 
Ao = 5 (im (GR5; dot-dashed). Again, the fiducial .model is included (solid), as are UVB and FIRB limits as in Figure El with FIRAS FIRB 
detection limits in dark, solid. See text and Table hi for details on models. 



